Local tensor network for strongly correlated projective states 



B. Beri and N. R. Cooper 
Theory of Condensed Matter Group, Cavendish Laboratory, J. J. Thomson Ave., Cambridge CBS OHE, UK 

(Dated: January 2011) 

The success of tensor network approaches in simulating strongly correlated quantum systems 
crucially depends on whether the many body states that are relevant for the problem can be encoded 
in a local tensor network. Despite numerous efforts, strongly correlated projective states, fractional 
quantum Ifall states in particular, have not yet found a local tensor network representation. Here 
we show that one can encode the calculation of averages of local operators in a Grassmann tensor 
network which is local. Our construction is explicit, and allows the use of physically motivated trial 
wavefunctions as starting points in tensor network variational calculations. 

PACS numbers: 73.43.-1,71. 27.+a 



Approximating ground states of strongly correlated 
quantum systems is one of the central challenges in con- 
densed matter physics. In recent years there have been 
several promising proposals aimed at tackling thisprob- 
lem for gapped spin[l|-y|, as well as fermionic[5|-|9| sys- 
tems with local interactions. Common to these ap- 
proaches is that they are all based on a class of states 
which can be parametrized with tensor networks. For a 
lattice system, the tensors live on the sites, which are con- 
nected by a set of links such that the indices of tensors at 
the ends of the links are contracted (possibly according 
to a metric defined by link tensors). An approximation 
of the ground state is obtained by treating the tensors as 
variational degrees of freedom. In order for tensor net- 
work approaches to be efficient, it is important that, for 
an accurate representation of the state, the number of 
degrees of freedom per tensor does not grow exponen- 
tially with the system size. Tensor networks with this 
property are local tensor networks. There exist various 
schemes for efficiently calculating averages of local oper- 
ators with respect to local tensor network states 
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Strongly correlated systems can realize gapped phases 
showing topological order, with features (ground state 
degeneracy, braiding statistics, etc.) that depend only 
on the topolo gy o f the configuration space, not its lo- 
cal details [13, Due to the topological nature of 
these phases, it is a nontrivial question whether one can 
capture topologically ordered ground states with local 
tensor networks. For spin systems, Refs. 
have shown that the ground states in a class of time- 
reversal invariant topological phases, so called string-net 
condensates 16|, ll7| , admit a local tensor network repre- 
sentation. The experimentally most accessible topologi- 
cally ordered states, the fractional quantum Hall (FQH) 
states, however, are fermionic, break time-reversal sym- 
metry. Their existing tensor network representations in 
the literature turn out to be nonlocal^, H, [11]. Does 
this mean that these states, and their spin-system de- 
scendants (such as the chiral spin liquid statejl^ Izoj ) 
lie outside of the scope of tensor network based varia- 
tional simulations? Fortunately, the results of this Let- 



ter demonstrate the contrary: we show that FQH states, 
or more generally, states from the so-called projective 
construction 2l|, [22] do admit a local tensor network rep- 
resentation. Tensor network based variational schemes 
should thus be able to find such ground states. 

To explain our idea, we start by first describing 
a non-interacting system: the many electron ground 
state l^*) of a gapped quadratic Hamiltonian H — 
1 V'/c)'i/ci('0i: where the indices k,l label 

sites of a lattice and •ipi^'^ are the fermionic field operators 
on these sites. (Writing iJ, in terms of a Bogoliubov de 
Gennes Hamiltonian h with a 2 x 2 electron-hole struc- 
ture for each kl is redundant for Hamiltonians without 
pairing, but it allows us to treat pairing states and states 
with fixed particle number on equal footing.) In the usual 
tensor network approaches, the tensor network represen- 
tation of 1 4") would be a preparatory step in express- 
ing averages of operators as traces over tensor networks. 
The key observation of our work is that, instead of first 
representing one can construct tensor networks for 
the averages directly. We consider operators of the form 



g^^(i/'])'^^ V'j^ , where kj,kj = 0, 1 and Ka is the 



set of sites j where A acts, i.e., where kj -\- kj 7^ 0. (A 
general local operator can be obtained as a linear com- 
bination of such products.) Using a result of Bravvi[23j| 
we can express (^'jAj^') as a Grassmann integral. 



X n (1) 

where Af is an A-independent normalization 
factor [i^l. (We omit Af henceforth, as it drops 
from (vI'|A|*)/(vI/|vl/).) Here 0^, 0j are Grass- 
mann variables, and iol^j^"^] is (j)j/\/2 / \/2) and 
ioltl^jipj] = exp{4>j(j)j)/2. The matrix h^^^ is the "flat 
band Hamiltonian" , which has the same eigenvectors as 
h, but has the eigenvalues Ej replaced by sgn{Ej)^ 



Using the Grassmann Leibniz rule (see e.g. Ref [2 
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Eq. ([T]) can be transformed from an integral over 
fermions (jjk living on lattice sites to one over fermions 
(j)K living on the legs (i.e., ends of links) K of sites k, 
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Here, the symbol Pq represents a projection of the re- 
sult of the integral to the term containing no Grassmann 
variables. The factors G^i — 9ki9ik with 



t — ift) 
gki cxp[--(0K,(/>i^) hli'i 



(3) 



belong to links kl between sites k and /. The precise form 
of the site factor Tj depends on the operators of A on j; 
in the case that j ^ Ka it is given by 



r, = (^d0^)(^d0^)-(/i(f)". 



(4) 



where the indices 11 refer to the electron- hole structure 
and J2 Jej denotes summation over the legs J of site j. 
Why did we go from Eq. ([T]) to Eq. ([2])? Because with 
Eq. (2)), (^'1^1^') takes the form of a tensor trace in an 
existing tensor network scheme, namely the scheme of 
Grassmann tensor networks introduced by Gu et al.^. 

Our construction has a saliently pleasant feature: it 
follows from h explicitly. This is to be compared with 
the usual situation in tensor network approaches, where 
if one has an explicit tensor network for a state j^") , often 
it is not derived from a Hamiltonian, but a Hamiltonian 
(for which j^*) is the ground state) is derived from it. 
The opposite direction is viable only numerically. The 
key feature of our construction, however, is rooted in 
the fact that for gapped systems, the matrix elements 
hlj decay exponentially as the function of |z — j|[25, ,28]. 
This implies that our tensor network is local. Indeed, in- 
troducing a cutoff length Z, the relative error introduced 
by neglecting a link longer than I is overestimated by 
exp(— Z/<^), where <^ is the decay length of h'-^\ The total 
number of links is N'^, thus iV^exp(— Z/^) largely over- 
estimates the error due the neglected links. Requiring 
that the error remains fixed implies that we have to scale 
the cutoff at most as I ^ InN. The number of legs at 
site j is estimated as l'^ ~ (In A'')'' in a d dimensional 
system. Since the number of variables M per leg is in- 
dependent of the system size, the number of degrees of 
freedom in Tj (equal to the number of coefficients in its 
expansion in d(j)K, etc. [9]) grows at most as 2[*^('"^) 1, 
i.e., quasi-polynomially. 

Solving quadratic Hamiltonians with tensor networks 
is certainly not a practical strategy. The usefulness of 
our construction lies elsewhere: as we now explain, it can 
be employed to describe an important subset of strongly 
correlated states, the states from the projective construc- 
tion. These states play an important role in several are- 
nas of strongly correlated quantum systems, such as the 



FQH effect, spin-liquids, high Tc superconductors, and 
such novel states of matter as the d = 3 fractional topo- 
logical insulator state [2§|. 

The projective construction defines strongly correlated 
trial ground states of the form 21 , 2^ 



(5) 



The state is written in the occupation number representa- 
tion, with Uj being the occupation of site j. The state \ Q) 
is the ground state of a gapped quadratic Hamiltonian h 
of P flavors of fermions (partons) living in a configura- 
tion space identical to the physical one. The annihilation 
operators of the partons are {cQj|a = 1 . . . P}. The state 
1 4') can be bosonic or fermionic. For bosonic (fermionic) 
states, the operators ijjj are even (odd) polynomials of 
the operators Caj. The operator ipj should not be con- 
fused with the field operator -0^ : the former acts in the 
parton Fock space, and is not a field operator, while the 
latter acts in the physical Fock space, and is a fermionic 
or bosonic field. 

A well known projective state is the Laughlin state at 
filling fraction = 1/3(33. Its wavefunction in position 
representation is 'i'{{zi}) ~ 4'i({zi})^, where "iiidzi}) 
is the wavefunction of the integer quantum Hall (IQH) 
state at 1/ = 1. In terms of Eq. ([5]), ^3 corresponds to 
P = 3, = Ci,C2jC3, and |e) = 1*1)2 \^i)3- This 

Laughlin state is a member of a class of states with|22| 



Qi— 1 



(6) 



where S* < P an integer. Depending on one can ar- 
rive at various FQH states, such as the Zk parafermion 
3E 



states 



31 



or the states based on the composite 
fermion picture [sj]- 

To use our tensor network construction for projective 
states, we need to convert averages of the form (*|yl|5') 
into averages with respect to the state \Q). (As before, 
we consider A 
we have 



rijGKyi (^1)'^^''^^^ ■) bosonic case, 



n,-!(n,- -t- kj - fc,)! 



\nj + kj - kj) (7) 



if Uj > kj and otherwise. This leads to (with kj = kj = 
for j i Ka) 

(vi/|A|M^) = (ein 



E :L mmiY 



19). 

(8) 

The key feature is that the expectation value of an op- 
erator factorizes into a product over different j-s. This 
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is essential for obtaining a local tensor network in the 
subsequent stages. 

In the fermionic case, the formula similar to Eq. ([7]) 
contains a sign factor depending on the occupations of 
sites j < I, preventing such a factorization. This can 
be remedied by using objects with anticommuting ingre- 
dients. More precisely, suppose one can find an opera- 
tor Rj which depends only on partons at j, and com- 
bines with ^jjj^ to mimic field operators in the following 
sense: {Rj,Rj>} — 0, and for / ^ j {_Rj, = and 
Rjijj]\0) = |0), i?j|0) = 0. One can then show that 



(vl/|^(^t^^)|vi,) = {Q\A{i>\R) 



n 



E 



19), (9) 



our construction can be built up straightforwardly. Just 
as in Eq. ([T]) , we have averages with respect to the ground 
state \Q) of a gapped quadratic fermion Hamiltonian. In 
both the boson and the fermion case, we have a prod- 
uct between {Q\ and |G) which we denote Ilj/i- The 
average will then have the same form as in Eq. ([T]), 

but with Uj^ifj) instead of Jlj-g/^^ ]■ The 

factor uj{fj) is obtained similarly as before: one takes 
fj{cl^j,Caj), expands it as a normal ordered polynomial 

in the parton operators. If in a term both cj^^- and Caj are 
present, they are brought next to each other. Single fac- 
tors Caj (Caj) are then substituted by ipaj/V^ {(j)aj/\/2), 
and the product c'ajCaj by exp[^„j0Qj]/2. 

The Grassmann Leibniz rule leads again to an expres- 
sion for (^I'l^l^') in terms of a tensor network of the form 
0. The link factors are now given by Gki =gkigik, 



where the ordering in the product is arbitrary, since the 
factors are even in parton operators. Since A{il}\ R) also 
factorizes, we now have an average over a product as we 
wanted. The order of factors matters only for j e Ka 
with odd kj + kj. The question is whether we can find 
such an Rj . We present here the solution of this represen- 
tation problem for the class of states described by Eq. ([6]) . 
In that case, Rj = ipj satisfies all the requirements, if x 
is chosen to satisfy X^^^^j Xai...asX*a,...as = '^f^-- 

Expressions (IH),® give a starting point from which 



gki=exp 



where /i^,^'' is a 2P x 2P block between partons at sites 
i and j of the fiat band Hamiltonian associated to \Q). 
The Grassmann variables 4>aK-:<PaK live on the leg of k 
ending the link Ik, etc. 

To obtain the site factors, one makes the expansion 



a/3 



, (10) 



'^(/j)exp 



Q/3 



(11) 



{TTa,j,TTa,j} 



where we introduced the expansion coefficients The product arranges the as in increasing order. A 

straightforward calculation leads to 



{Traj-.TTaj} 



{7r„j,7r„j} 



n (1 - 2^„,)'"*°^(E^'^"^)'""°'(E^^"^)'~*° 



(12) 



This tensor network for the states of the projective con- 
struction has all the discussed advantages: it is formu- 
lated in the framework of an existing scheme, it is ex- 
plicitly constructed, and most importantly, it is local. 
Eqs. ([2]), ([To)) and p2| are the main results of the Let- 
ter. 



We now use our results to present the link and site 
factors for some concrete projective states. We first 
consider a fermionic FQH state I^'cf) on the lattice 
from the composite fermion series with Hall conductance 
CTi2 =m/(l + 2m) (with m integer). It is a generalization 
of the v = 1/3 state, obtained by replacin g on e of the 
|*i)-s by I^I'm), an IQH state with ai2 = Trtjsi. IsHj. Its 



tensor network representation has 



(13) 

where hm is the fiat band Hamiltonian obtained from 
the IQH Hamiltonian hm realizing the state |^'m). The 
site factors for j ^ Ka are given by 

^i^7(n°^aj+X]'^ai)' d(j)a.jd4>aK~H°'j- (14) 



If the operator content of A on site j is ipj , we have 

^ Q ,7ei 



(15) 
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The other types of site factors can be found similarly. 

In the CF state, tpj was a product of c-s. As pointed 
out in Ref. Q , calculating averages with respect to states 
with such i/jj is also possible by variational Monte Carlo 
techniques, but these techniques fail once 'ijjj is a sum of 
products of c-s. It is thus important that our construc- 
tion includes the latter case as well. As an illustration we 
consider (the lattice version of) the bosonic Pfaffian state. 
In terms ofEq. (0, this state has |e) = j^-i) j^-i) l^fi) 
and "ipj — cijCij — C3jC2j^2^. Its tensor network represen- 
tation has 



gki = n 



(fl)Nll 



(16) 



The site factors, for example for j (f. Ka, are given by 
= ]^[n("ai + 1) + 4(aija4j + I)(a2-,a3j -|- 1) 

a 

-4 ^ {d(f>ijd(l)2Kd(f>3Ld(j>iM+d(f>ijd(j>2Kd(f>3Ld(l)4M)]- 
J,K,L,Mej 

Once the site and link factors Tj and Gki are known, 
the calculation of (5'|A|^') amounts to the evaluation of 
the tensor trace Eq. ([2]). In tensor network approaches, 
evaluating tensor traces is a challenge by itself, which re- 
quires approximations. Depending on the concrete ten- 
sor network framework at hand, one can choose from 
existing approximation schemes to tackle this task[l|-Q. 
As Eq. ([2]) is formulated entirely in terms of the frame- 
work of Ref. i, it can be subjected to the approximation 
scheme developed there, the Grassmann tensor entangle- 
ment rcnormalization group. The concrete implementa- 
tion of this rcnormalization procedure for our networks 
is however beyond the scope of this Letter. 

In conclusion, we have demonstrated that strongly cor- 
related projective states can be represented using local 
Grassmann tensor networks. The results are valid for the 
whole scope of the projective construction, which include 
FQH states, spin-liquids, fractional topological insula- 
tors, etc. On the level of principles, our work shows that 
strongly correlated quantum systems where the ground 
state is close to a projective state can be meaningfully ap- 
proached using tensor network based algorithms. From a 
computational point of view, our development brings the 
prospects of variational simulations with projective trial 
states qualitatively closer: first, it breaks an exponen- 
tial bound faced in previous efforts in representing such 
states via tensor networks; and, second, it does so with- 
out giving up on the explicit nature of the representation. 
This allows to use physically motivated trial wavefunc- 
tions as starting points, which could speed up variational 
calculations dramatically. 

This work was supported by EPSRC Grant 
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